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ABSTRACT 

Using an efficient numerical scheme that exploits spatial symmetries and spin parity, we have 
obtained the exact low-lying eigenstates of exchange Hamiltonians for the high nuclearity spin 
clusters, Mni2, Fes and V15. The largest calculation involves the Mni2 cluster which spans a Fock 
space of a hundred million. Our results show that the earlier estimates of the exchange constants 
need to be revised for the Mni2 cluster to explain the level ordering of low-lying eigenstates. In 
the case of the Fes cluster, correct level ordering can be obtained which is consistent with the 
exchange constants for the already known clusters with butterfly structure. In the V15 cluster, we 
obtain an effective Hamiltonian that reproduces exactly, the eight low-lying eigenvalues of the full 
Hamiltonian. 



PACS numbers: 75.50.Xx, 61.46.-l-w 



I. INTRODUCTION 

The synthesis of high nuclearity transition metal complexes has provided a new dimension to 
the field of nanomagnetism.^ Many interesting phenomena have been observed in these systems. 
Amongst the most exciting are the observation of quantum resonance tunneling and quantum 
interference^ in some of these clusters. For example, in the case of the Mni2 cluster, the ground 
state with total spin Sg = 10 of the exchange Hamiltonian, under the influence of a large single 
ion anisotropy gives rise to a manifold of doubly degenerate states with nonzero Ms values, with 
Ms = ±10 being the lowest energy states. The application of a magnetic field splits the degeneracy 
of the Ms = ±10 states. Varying the magnetic field brings states with |Afs| ^ 10 closer in energy 
to the higher of the two states with \Ms \ — 10. The weak spin dipolar interactions that exist in the 
system can connect these nearly degenerate states with different Ms values, leading to tunneling 
between the states. This is reflected in experiments as jumps in magnetization in the magnetization 
vs magnetic field plots, whenever the resonance condition is satisfied and as plateaus for off- 
resonance field values. Similar plateaus are also observed in the Vi^ cluster'^, although the reason 
for the plateaus in this system is qualitatively different. The quantum interference phenomena 
observed in the Fes cluster is because the paths connecting the Ms = +10 and Ms = —10 could 
interfere in the presence of a magnetic field, leading to an oscillation in the tunneling probabilities.'' 

These clusters, at a very basic level are characterized by multidentate ligands interconnecting the 
transition metal ions. In the clusters, a given magnetic ion has exchange interactions of either sign 
with several of its neighbors. Thus, these magnetic clusters often correspond to spin frustrated 
systems. Because of the rather complex exchange pathways which exist in these systems, it is 
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difBcult to predict a priori even the sign of the exchange constant, let alone its magnitude.^ Since 
the site symmetry at the magnetic ions is also usually low because of the multidentate ligands 
in the system, the orbital degeneracy would be lifted leading to weaker ferromagnetic exchange 
interactions. Thus one should expect low-spin ground states in these systems. However, because 
of the frustrations in the exchange pathways, even weak ferromagnetic interactions could lead to 
higher spin ground states, albeit with rather low spin excitation gaps. It is indeed interesting 
that the Mni2 cluster has ferrimagnetic ordering^ of the spins in its ground state for this reason. 
The high-spin ground state observed in Feg, is, however, attributable entirely to frustration in the 
antiferromagnetic interactions.^ 

In the case of the Mni2 cluster, while the ground state spin as well as the lowest excitation gap 
is established experimentally, it is not at all clear what the magnitude and sign of the exchange 
interactions in the cluster are. In an earlier study^, in order to simplify the calculations, each 
strongly coupled Mn^^^ — Mn^^ pair was replaced by a composite spin-1/2 object. The resulting 
model was studied for three different sets of exchange constants. It was observed that the ordering 
of the energy levels is very sensitive to the variations in the exchange constants. 

In the case of the Fcg cluster, while model exact calculations® were possible because of the smaller 
dimensions of the Hilbert spaces, the exchange parameters used were very different from those that 
have been determined recently.® Considering the sensitivity of the ordering of the energy levels to 
values of the exchange constants, it is desirable to redo the calculations using revised estimates of 
the exchange constants. 

The simplest cluster that can be studied exactly is the V15 cluster. In this cluster it is found 
that eight low-lyings states are well separated from the rest of the spectrum.^ Most of the low- 
temperature properties are determined by these eight low-lying states. To undertake serious study 
of the magnetization behavior of the cluster under the influence of an applied magnetic field, it is 
necessary to construct an accurate model Hamiltonian for these states. 

In this paper, we have studied the low-lying states of Mni2, Fes and V15 clusters using exact 
diagonalization of the corresponding exchange Hamiltonians. We show that in the case of Mni2, 
earlier estimates of the exchanges constants fail to provide the ground and excited state spin 
quantum numbers in agreement with experiments. We have also estimated the most likely exchange 
constants that give good agreement with experiments. In the case of the Fes cluster, we have 
studied the properties of the cluster for more recent estimates of the exchange constant. In the 
case of V15 we have obtained an effective Hamiltonian for the low- lying states that reproduces the 
energy level ordering of the eight low- lying states exactly. Such a Hamiltonian would be important 
in the context of hysteresis studies of the system at low temperatures. 

In the next section, we outline the numerical method for obtaining the low-lying states of ex- 
change coupled spin systems which span large Hilbert spaces. In Sec. 3, we discuss the results 
obtained from this technique for the three magnetic clusters mentioned above, and we summarize 
our results in Sec. 4. 



II. MODEL HAMILTONIAN AND COMPUTATION DETAILS 



The model Hamiltonian employed in these studies is the isotropic exchange Hamiltonian involv- 
ing exchange interactions between nearest neighbors, 

H = E Jij Si ■ Sj , (1) 

<ij> 
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where the exchange interaction takes the values dictated by experimental studies of structure 
and magnetic properties. The total dimensionality of the Fock space of the cluster is given by 

n 

Dp = n {2Si + 1) , (2) 
1=1 

where n is the total number of spins in the cluster and Si is the spin on each ion. In the case of the 
Mni2 cluster consisting of eight spin-2 ions and four spin-3/2 ions, the Fock space dimensionality is 
a hundred million. Specializing to a given total Ms leads to Hilbert space dimensionalities which 
are lower than the Fock space dimensionality. In the case of the Mni2 cluster the Ms = space 
has a dimensionality of over eight million (8,581,300). The major challenge in exact computation 
of the eigenvalues, and properties of these spin clusters lies in handling such large basis and the 
associated matrices. While the dimensions look overwhelming, the matrices that represent the 
operators in these spaces are rather sparse. Usually, the number of nonzero elements in a row is of 
the order of the number of exchange constants in the Hamiltonian. This sparseness of the matrices 
allows one to handle fairly large systems. However, in the case of spin problems, generating the 
basis states and using the symmetries of the problem is nontrivial. 

The isotropic exchange Hamiltonians conserve the total spin, S, besides the 2-component of 
the total spin, Ms- Besides these symmetries, the geometry of the cluster also leads to spatial 
symmetries which can often be exploited. The simplest way of generating basis functions which 
conserve total spin is the VB method that employs the Rumer-Pauling rulc.^° It is quite easy to 
generalize the Rumer-Pauling rules to a cluster consisting of objects with different spins to obtain 
states with desired total spin, 5*. However, setting up the Hamiltonian matrix in such a basis can be 
computer intensive since the exchange operators operating on a " legal" VB diagram (diagram that 
obeys Rumer-Pauling rules) can lead to "illegal" VB diagrams, and resolving these "illegal" VB 
diagrams into "legal" diagrams would present the major bottle-neck. Indeed, the same difficulty is 
encountered when spatial symmetry operators operate on a VB function. Thus, the extended VB 
methods are not favored whenever one wishes to apply it to a motley collection of spins or when 
one wishes to exploit some general spatial symmetries that may exist in the cluster. 

Usually, in frustrated spin systems, it is important to partition the spaces into different total 
spin spaces because of the usually small energy gaps between total spin states which differ in S by 
unity. To avoid the difficulties involved in working with total spin eigenfunctions, we exploit parity 
symmetry in the systems. The parity operation involves changing the z-component of all the spins 
in the cluster from Mg^ to —Ms^. There is an associated phase factor with this operation given by 
s s- 

(— 1) '"'"""Zji *. The isotropic exchange operator remains invariant under this operation. If this 
symmetry is employed in the Ms = subspace, the subspace is divided into "even" and "odd" 
parity spaces depending upon the sign of the character under the irreducible representation of the 
parity group. The space which corresponds to even (odd) total spin we call the even (odd) parity 
space. Thus, employing parity allows partial spin symmetry adaptation which separates successive 
total spin spaces, without introducing the complications encountered in the VB basis. However, 
the VB method can lead to complete factorization of the spin space leading to smaller complete 
subspaces. 

In the Mni2 cluster, besides spin symmetries, there also exists spatial symmetries. The topology 
of the exchange interaction leads to a C4 point group symmetry. At first sight, this point group 
appears to present difficulties because the characters in the irreducible representation are in some 
cases complex. This could lead to complex basis functions. This, however, can be avoided by 
recognizing that in the C„ group, states with wave vectors k and —k are degenerate in the absence 
of an external magnetic field. We can therefore construct a linear combination of the k and —k 
states which is real. The symmetry representations in the C4 group would then correspond to 
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the labels A, B and E, with the characters in the E representation given by 2cos(rfc) under the 
symmetry operation CJ, with k = 7r/2. The parity operation commutes with the spatial symmetry 
operations, and the full point group of the system would then correspond to the direct product of 
the two groups. Since both parity and spatial symmetries can be easily incorporated in a constant 
Ms basis, we do not encounter the difficulties endemic to the VB theory. 

The generation of the complete basis in a given Hilbcrt space requires a simple representation 
of a state on the computer. This is achieved by associating with every state a unique integer. 
In this integer, we associate rij bits with spin Si, such that rii is the smallest integer for which 
2"' > 2si + 1. In the integer that represents the state of the cluster, we ensure that these rii bits 
do not take values which lead to the nj-bit integer value exceeding 2si + 1. For each of the allowed 
bit states of the rij-bit integer, we associate an M5, value between — Sj and Sj. For a spin cluster 

n 

of n spins, we scan all integers of bit length N = T, rii and veriiy if it represents a basis state with 

i=l 

the desired Ms value. In Fig. 1, we show a few basis functions with specified Ms value for some 
typical clusters along with their bit representations and the corresponding integers. Generation of 
the basis states is usually a very fast step, computationally. Generating the basis as an ordered 
sequence of integers that represent them also allows for a rapid generation of the Hamiltonian 
matrix elements as will be seen later. 

Symmetrization of the basis by incorporating parity and spatial symmetries involves operating 

on the constant Ms basis by the symmetry operators. Since spatial symmetry operators exchange 
the positions of equivalent spins, every spatial symmetry operator operating on a basis function 
generates another basis function. Every symmetry operator can be represented by a correspondence 
vector whose i^^ entry gives the state that results from operating on the i**^ state by the chosen 
operator. This is also true for the parity operator, in the Ms = subspace. The symmetry 
combinations can now be obtained by operating on each state by the group theoretic projection 
operator, 

^r. = l^XvmR (3) 

on each of the basis states. Here Fj is the i^^ irreducible representation, R is the symmetry 
operation of the group, h is the order of the symmetry group, and Xri(^) is the character under R 
in the irreducible representation F^. The resulting symmetrized basis is overcompletc. The linear 
dependencies can be eliminated by a Gram-Schmidt orthonormalization procedure. However, in 
most cases, ensuring that a given basis function does not appear more than once in a symmetrized 
basis is sufficient to guarantee linear independence and weed out the linearly dependent states. 
A good check on the procedure is to ensure that the dimensionality of the symmetrized space 
agrees with that calculated from the traces of the reducible representation obtained from the 
matrices corresponding to the symmetry operators. Besides, the sum of the dimensionalities of the 
symmetrized spaces should correspond to the dimensionality of the unsymmetrized Hilbert space. 

The generation of the Hamiltonian matrix is rather straightforward and involves operation of 

the Hamiltonian operator on the symmetry adapted basis. This results in the matrix SH, where 
S is the symmetrization matrix representing the operator P^^ and H is the matrix whose elements 
hij are defined by 

H\i>=^hij\j> . (4) 
3 

The states i correspond to the unsymmetrized basis functions. The Hamiltonian matrix in the 
symmetrized basis is obtained by right multiplying the matrix SH by S^. The symmetric Hamilto- 



4 



nian matrix is stored in the sparse matrix form and the matrix eigenvalue problem is solved using 
the Davidson algorithm. 

Computation of the properties is easily done by transforming the eigenstate in the symmetrized 
basis into that in the unsymmetrized basis. Since the operation by any combination of spin 
operators on the unsymmetrized basis can be carried out, all relevant static properties in different 
eigenstates can be obtained quite simply. 



III. RESULTS AND DISCUSSION 

We have solved the exchange Hamiltonian exactly for the Mni2, Fes and V15 clusters using 
the above mentioned method. We have obtained the eigenvalues and various properties of the 
eigenstates such as spin densities and spin-spin correlation functions for these clusters. In what 
follows, we will discuss these in detail. 



A. Mni2Ac Cluster 

In Fig. 2 we show the geometry and the exchange parameters for this cluster. The crystal 

structure suggests that the exchange constant Ji is largest and antifcrromagnctic in nature. -'^'^ Based 
on magnetic measurements, it has been suggested that Ji has a magnitude of 215K. The other 
magnitude and sign of the other exchange constants are based on comparisons with manganese 
systems in smaller clusters. '^-'^ It has been suggested that the exchange constant J2 and J3 are 
antifcrromagnctic and have a magnitude of about 85K. However, for the exchange constant J4, 
there is no concrete; estimate, either of the sign or of the magnitude. In an earlier study, the Mn-'^''^-'^ 
- Mn^^ pair with the strongest antiferromagnetic exchange constant was replaced by a composite 
spin-1/2 object'', and the exchange Hamiltonian of the cluster was solved for three different sets of 
parameters. It was found that the ordering of the energy levels were very sensitive to the relative 
strengths of the exchange constants. In these studies, J4 was set to zero and the low-lying excited 
states were computed. Besides, only states with spin S up to 10 could be obtained because of the 
replacement of the higher spin ions by composite spin-1/2 objects. 

In our calculation, we have dealt with all the magnetic ions in the cluster and using symmetry, 
we have factored the Ms = Hilbert space into the six symmetry subspaces. The dimensionalities 

of the different subspaces is given in Table I. Wc have obtained low-lying eigenstates in each of 
these sectors and determined the total spin of the state by explicitly computing the expectation 
value of the S"^ operator in the state. 

Our results for the low-lying states are shown in Table II. We note that none of the three sets 
of parameters studied using an effective Hamiltonian, gives the correct ground and excited states, 
when an exact calculation is performed. It appears that setting the exchange constant J4 to zero, 
cannot yield an S = 10 ground state (Table II, cases A, B and C). When J3 is equal to or slightly 
larger than J2 (cases A and B, Table II), we find a singlet ground state, unlike the result of the 
effective Hamiltonian in which the ground state has 5 = 8 and 5 = respectively. The ground 
state has spin 5 = 6, when J3 is slightly smaller than J2 (case C, Table II). In all these cases, the 
first few low-lying states are found to lie within 20K of the ground state. 

When we use the parameters suggested by Chudnovsky^^ (case D, Table II), we obtain an 5 = 10 
ground state separated from an 5 = 9 first excited state by 223K. This is followed by another 5 = 9 
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excited state at 421K. Only when the exchange constant J4 is sufficiently strongly ferromagnetic 
(case E, Table II), do we find an 5 = 10 ground state with an 5 = 9 excited state separated from 

it by a gap of 35K, which is close to the value inferred indirectly from experimental results. The 
second higher excited state has 5 = 8, and is separated from the ground state by 62K. 

We have explored the parameter space a little ftirther by varying J3 and J4, to see the eflFect of 
these exchange constants on the ordering of the energy levels. We find that for | J3I = | J4I and J3 
antiferromagnetic but J4 ferromagnetic, the ground state is always 5 = 10 (Table 111, cases C, D 
and E); the first and second excited states are 5 = 9 and 5 = 8, respectively. The lowest excitation 
gap decreases slowly with increasing magnitude of the exchange constants. 

We find that the spin of the ground state is very sensitive to J4, for fixed value of J3. In 
the case where J4 weakly ferromagnetic (Table III, case B), we obtain an 5 = ground state, 
J4 weakly antiferromagnetic we obtain an 5 = 4 ground state (Table III, case A). This shows 
that frustrations play a dominant role. If J3 is also made ferromagnetic, the role of frustration is 
considerably reduced. 

In Figs. 3(a) and 3(b), we show the spin density-^^ for the Mni2 cluster in the ground state for 
the 5 = 10, Ms = 10 state. While the manganese ions connected by the strong antiferromagnetic 
exchange show opposite spin densities, it is worth noting that the total spin density on these two 
ions is 0.691, well away from the value of 0.5 expected, if these ions were indeed to form a spin-1/2 
object. We also note that the spin density at the manganese ion in the middle of the crown is much 
larger than that at the corners. The spin density in the excited state 5 = 9, Ms = 9, also has a 
similar distribution, although in this state, the symmetry of the spin Hamiltonian is apparently 
broken [Fig. 3(b)]. The corner ions in the crown no longer have the same spin densities; one pair 
of opposite corner ions have larger spin densities than the other corner pair. However, since this 
state is doiibly degenerate, there is another state in which the spin densities are related to the spin 
densities of this state by a 90" rotation. In any experiment involving this state, only an arbitrary 
linear combination of the two spin densities will be observed. Note also that the large differences 
in the spin densities for the closely lying excited states is an indication of the large degree of spin 
frustration in the system. 

The small energy gap (35K) between the 5 = 10 ground state and the 5 = 9 lowest excited state 

seems to indicate that, if the g factors of the Mn ions in the core and crown are different then, 
an applied magnetic field should mix these spin states. Such a mixing would then be reflected in 
the quantum resonance tunneling experiments. However, it appears that the experiments are well 
described by the 5 = 10 state alone. This is what we should expect from the symmetry of the two 
low-lying states. We note that the ground state has A symmetry while the lowest excited state 
has a E symmetry. These two states cannot be mixed by any perturbation that retains the C4 
symmetry of the cluster. 



B. Fes Cluster 

The Fes cluster is shown in Fig. 4. Each of the Fe ions has a spin of 5/2 and the ground state of 
the system has a total spin 5 = 10, with 5 = 9 excited state separated from it by about 20K. All 
the exchange interactions in this system are expected to be antiferromagnetic. While the structure 
of the complex dictates that the exchange interaction J2 along the back of the butterfly should be 
considerably smaller than the interaction Ji across the wing-'^^, in earlier studies it was reported 
that such a choice of interaction parameters would not provide a 5 = 10 ground state. ^ 

We have carried out exact calculations of the eigenstates of the Feg cluster using four sets of 
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parameters; the last set of parameters (Case D) are taken from Ref. [17]. In cases A, C and D, 
J2 is very much smaller than Ji. We find that in all the four cases, the ground state has a spin 
S' = 10 and the lowest excited state has spin S = 9. One of the main differences we find amongst 
the four sets of parameters is in the energy gap to the lowest excited state (Table IV) . For the set 
of parameters used in the earlier study, this gap is the lowest at 3.4K (case B). For the parameter 
sets A, C and D®, this gap is respectively 13. IK, 39. 6K and 42. 4K. While in cases A, C and D, the 
second excited state has spin-9, in case B, this state has spin-8. 

The spin densities in all the four cases for both the ground and the excited state are shown in 

Figs. 5(a) to 5(h). The spin densities in all cases arc positive at the corners. In cases A and B, 
the spin density is negative on the Fe ions on the backbone, and is positive on the remaining two 
Fc sites. ^^'-^^ However, in cases C and D, the negative and positive spin density sites for the Fe 
ions in the middle of the edges are interchanged. This is perhaps due to the fact that in cases A 
and B, the exchange constant J3 is less than J4, while in cases C and D, this is reversed. Thus, 
a spin density measurement can provide relative strengths of these two exchange constants. In all 
the cases, the difference between the spin densities in the ground and excited states is that the 
decrease in the spin density in the excited state is mainly confined to the corner Fe sites. Note 
that the spin densities in cases C and D are almost the same, although the excitation gaps are 
significantly different and in proportion to the differences in the exchange constants. 

We should note here that the spin densities presented by us are expectation values of the site 
operators in the S = 10, Ms = 10 and 5* = 9, M5 = 9 states. However, experimental values are 
obtained not only at the Fe sites but also at the spin polarized neighboring ligand atom sites. We 
therefore use the experimental results only as a guideline for locating the negative spin density Fe 
sites which are sensitive to the set of exchange constants used in the calculation. 

The Peg cluster is quite different from the Mni2 cluster in the following sense. In the Fcs cluster, 
we have excited states of spin S = 8 and S = 9 which have the same symmetry as that of the 
S = 10 ground state. Furthermore, the total splitting of the ground state due to the anisotropic 
terms arising in the system due to spin-dipolar interactions is larger than the energy gaps with the 
5 = 8 and 5 = 9 states of the same spatial symmetry as the ground state. Thus, if the g factors of 
the Fe ions on the backbone of the butterfly are different from those on the wings, then an applied 
magnetic field could lead to mixing between the different spin states. We expect this to provide 
an additional mechanism for quantum resonance tunneling in the Feg cluster. 



C. Vi5 Cluster 

The simplest cluster to study is the V15 cluster, since each of the ions has a spin of half. The 
interesting aspect of the V15 cluster is that the three spins sandwiched between the hexagons [Fig. 
6] have no direct spin-spin interactions. All the interactions shown in Fig. 6 are antiferromagnetic 
and the spin system is frustrated. The eigenstates of this system consist of eight states correspond- 
ing to the triangle spins (i.e., three spin-1/2 sites) which are split off from the rest of the spectrum. 
A combination of three two-spin interactions which retains the C3 symmetry of the molecule is 
sufficient to account for such a spectrum. We find that the effective Hamiltonian is given by 

H,p_,p = e / + a (5i • 52 + 52 • 53 + 53 • 5i) , (5) 

where e = —4.78187 and a = 0.02015 in units of the exchange Ji. This Hamiltonian reproduces 
the eight low-lying eigenstates of the full exchange Hamiltonian to numerical accuracy. 
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The spin density distribution in one of the S = 1/2, Ms = 1/2 ground states as weh as in the 
S = 3/2, Ms = 3/2 excited state is shown in Figs. 7(a) and 7(b). We find that the spin densities 
on the hexagons are negHgible. The total spin density in the triangle for the 5=1/2 state is 
nearly equal in value to that a free electron spin. The S = 3/2, Ms = 3/2 state also has almost 
equal spin densities at all three sites of the triangle, nearly equal to that of free spins. These 
observations suggest that describing the low-energy spectrum of this system by the triangle spins 
is quite appropriate. 



IV. SUMMARY 



To conclude, using a bit representation of the spin states of a spin cluster combined with ex- 
ploitation of spatial symmetry and spin parity, we are able to obtain model exact solutions for 
exchange Hamiltonians whose Fock space spans up to a hundred million states. Our studies on 
the Mni2Ac cluster cluster shows that the earlier effective Hamiltonian studies wrongly estimated 
the exchange constants. The new exchange constants give the correct spin for the ground state as 
well as the correct ordering of the low-lying excited states. The spin densities in the cluster also 
support the fact that the effective Hamiltonian in earlier studies does not accurately represent the 
cluster. The studies on Feg cluster shows that the correct ground state as well as the first excited 
state can be obtained by using a set of exchange constants consistent with the butterfly structure 
known in related systems, unlike what was concluded in an earlier study. The interesting feature 
of the Fcs spin densities is that the sites which have negative spin densities depend on the relative 
strengths of the exchange constants on the sides of the butterfly. In the case of the V15 cluster, 
we flnd that the eight low-lying energy levels can be fitted to an effective Hamiltonian describing 
three spin-1/2 sites. 
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Table I: Dimensions of different symmetry subspaces in Ms=0 sector for Mni2Ac. 



Symmetry 


Dimension 




1074087 


°A 


1071537 




1074037 


°B 


1071587 




2142526 




2147526 



Table II: Low-Lying states of Mni2Ac, relative to the ground state for the parameters in question. 
Entries in parenthesis in cases A, B and C correspond to the effective Hamiltonian results of Scssoli 
et alJ Case D corresponds to the parameters suggested by Chudnovsky.^^ The parameters corre- 
sponding to different cases are: case (A) Ji=225K, J2=90K, J3=90K, J4=0K; case (B) Ji=225K, 
J2=90K, J3=93.8K, J4=0K; case (C) Ji=225K, J2=90K, J3=86.2K, J4=0K; case (D) Ji=215K, 
J2=85K, J3=-85K, J4=-45K; case (E) Ji=215K, J2=85K, J3=85K, J4=-64.5K. AU the energies 
are in K. 



Case A 


Case B 


Case C 


Case D 


Case E 


State S E(K) 


State S E(K) 


State S E(K) 


State S E(K) 


State S E(K) 


0.0 


0.0 


6 0.0 


""A 10 0.0 


M 10 0.0 


(8) 


(0) 


(10) 






"E 1 10.8 


"E 1 16.2 


"E 1 15.5 


"E 9 223 


"E 9 35.1 


(9) (6.4) 

(10) (6.4) 


(8) (1.4) 


(8) (2.7) 






°B 1 19.8 


°B 1 20.0 


"B 1 19.6 


°B 9 421.2 


«B 8 62.1 


(0) (6.8) 




(9) (5.0) 






"A 2 24.7 


M 2 30.5 


M 2 23.8 


°B 9 425.1 


"E 7 82.4 


"E 3 39.0 


"B 4 58.4 


°E 1 28.8 


"B 8 439.5 


"A 6 99.7 


"E 2 49.9 


2 60.9 


''B 6 53.6 


"B 8 443.7 


"B 102.0 


4 57.1 


"A 3 64.3 


"B 6 54.4 


8 458.1 


M 2 121.0 


^B 8 57.8 


^E 2 80.0 


^B 8 57.2 


°A 11 573.4 


"B 1 133.3 


''B 2 57.8 


°A 3 88.1 


''E 2 63.0 


"E 9 583.8 


^B 2 177.1 


"B 3 78.4 


M 6 88.3 


"A 3 77.0 


"E 8 632.8 


M 3 211.3 


°B 3 86.8 


"B 3 112.8 


"B 3 85.3 


°A 9 640.5 


°A 3 220.8 


"A 6 105.7 


°B 5 114.6 


•"E 2 86.1 


"E 8 658.3 


4 249.9 


"B 3 113.4 


"B 5 V>6A 


' .4 6 97.1 


' A 8 767.1 


"B 5 278.5 


"E 4 117.3 


°A 1 165.2 


''A 6 98.2 


"B 8 807.6 


M 7 332.1 


"B 5 154.2 


"A 1 181.6 


"B 3 112.2 


M 8 815.8 


M 7 340.8 
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Table III: Low Lying states of Mni2Ac. The parameters corresponding to different cases are: case 
(A) Ji=215K, J2=85K, J3=85K, J4=45K; case (B) Ji=215K, J2=85K, J3=85K, J4=-45K; case 
(C) Ji=215K, J2=85K, J3=64.5 K, J4=-64.5K; case (D) Ji=215K, J2=85K, J3=85K, J4=-85K; 
case (E) Ji=215K, J2=85K, J3=45K, J4=-45K. All the energies are in K. 



Case A 


Case B 


Case C 


Case D 


Case E 


State S E(K) 


State S E(K) 


State S E(K) 


State S E(K) 


State S E(K) 


' B 4 ().() 


' B U 0.0 


' .1 10 0.0 


' .1 10 0.0 


' .1 10 0.0 


"A 4 9.1 


°E 1 12.3 


°E 9 73.7 


"E 9 67.7 


°E 9 80.1 


3 9.4 


M 2 22.9 


"B 8 135.1 


^B 8 121.2 


^B 8 149.8 


^B 4 18.2 


•^B 1 27.6 


"E 7 186.1 


"E 7 165.2 


M 8 191.0 


"A 2 32.4 


"E 3 28.9 


M 8 196.0 


M 6 201.2 


"E 7 210.0 


"B 5 49.4 


^B 4 34.1 


''A 6 227.8 


M 8 206.5 


M 6 260.0 


^A 6 50.0 


M 10 36.5 


^B 4 283.5 


^B 4 247.7 


"B 4 329.8 


4 55.4 


8 37.8 


°B 1 323.0 


% 1 282.5 


"B 9 346.8 


"A 3 68.2 


<'E 2 67.2 


"E 2 364.0 


2 330.2 


°B 9 370.7 


"A 3 70.2 


°A 3 100.1 


°A 3 391.8 


M 3 365.2 


°B 1 515.8 


"B 3 71.4 


°A 3 119.5 


"A 3 401.6 


°A 3 375.0 


''E 8 400.3 


°A 3 76.6 


''A 4 140.0 


''E 4 420.6 


'^E' 4 401.9 


^E 2 413.8 


"B 2 255.2 


°B 3 161.8 


°B 9 426.3 


°B 11 421.0 


M 5 424.2 


«B 2 257.2 


°B 5 172.8 


"B 5 434.9 


°B 5 425.5 


M 3 432.5 



Table IV: Energies (in units of K) of a few low-lying states in Fes. The exchange constants 
corresponding to the various cases are: case (A) Ji = 150K, J2 = 25K, J3 = 30K, J4 = 50K; case 
(B) Ji = 180K, J2 = 153K, J3 = 22.5K, J4 = 52.5K; case (C) Ji = 195K, J2 = 30K, J3 = 52.5K, 
J4 = 22.5K; case (D) Ji = 201K, J2 = 36.2K, J3 = 58.3K, J4 = 26.1K. 



Case A 


Case B 


Case C 


Case D 


State S E(K) 


State S E(K) 


State S E(K) 


State S E(K) 


"A 10 0.0 
°B 9 13.1 
M 9 26.1 
M 8 27.3 


"A 10 0.0 
"B 9 3.4 
M 8 10.2 
°B 7 20.1 


"A 10 0.0 
"A 9 39.6 
"B 9 54.2 
°B 9 62.4 


''A 10 0.0 
"A 9 42.4 
°B 9 58.8 
°B 9 69.4 
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Figure Captions 



1. Representative Mg = state in (a) 6 spin-1/2 cluster, (b) Mni2Ac cluster with first four sites 

each having spin S = 3/2 and the remaining eight sites each having spin S = 2. Numbers in 
parenthesis correspond to the Mg = value at the site. The bit representations as well as the 
integer values are given just below the diagrams. 

2. A schematic diagram of the exchange interactions between the Mn ions in the Mni2Ac molecule. 

3. Spin density of Mni2Ac for parameter values: Ji= 215K, J2=85K, J3=85K and J4=-64.5K. (a) 
Spin density for ground state (S=10,Ms=10). (b) Spin density for 1*** excited state (S=9, Ms=9). 

4. A schematic diagram of the exchange interactions between the Fe ions in the Feg molecule. 

5. Spin density of Feg for parameter values: Ji = 150K, J2 = 25K, J3 = 30K, J4 = 50K. (a) Spin 
density for ground state (S=10,Ms=10). (b) Spin density for 1** excited state (S=9, Ms=9). 
Spin density of Fes for Ji = 180K, J2 = 153K, J3 = 22. 5K, J4 = 52. 5K parameter values, (c) Spin 
density for ground state (S=10,Ms=-10). (d) Spin density for 1^* excited state (S=9,Ms=9). 
Spin density of Fes for Ji = 195K, J2 = 30K, J3 = 52. 5K, J4 = 22. 5K parameter values, (e) Spin 
density for ground state (S=10, Ms=10). (f) Spin density for 1** excited state (S=9,Ms =9). Spin 
density of Fes for Ji = 201K, J2 = 36. 2K, J3 = 58. 3K, J4 = 26. IK parameter values, (g) Spin 
density for ground state (S=10, Ms=10). (h) Spin density for 1** excited state (S=9,Ms =9). 

6. A schematic diagram of the exchange interactions between the V ions in the V15 molecule. 

7. Spin density of V15 for Ji = 800K, J2 = 300K, J3 = 150K parameter values, (a) Spin density 
in one of the ground states (S=0.5,Ms=0.5). (b) Spin density for excited state (S=1.5, Ms= 1.5). 
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Fig. 1 
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Fig. 4 
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Fig. 5(b) 
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Fig. 5(c) 
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Fig. 5(d) 
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Fig. 5(e) 
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Fig. 5(f) 
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Fig. 5(h) 
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Fig. 7(a) 
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Fig. 7(b) 
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